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1. Introduction 



Over the last decade, event shape distributions in e + e~ annihilation have provided 
some of the most precise tests of QCD. Most of the attention has focussed on three-jet 
event shape variables such as thrust, wide jet broadening, heavy hemisphere mass, C 
parameter etc.. These variables may be generated by single gluon emission from the 
underlying qq pair in the e + e~ annihilation process which first contributes at 0(a s ). 
The analysis of such events is rather sophisticated. In the first instance, fixed order 
QCD perturbation theory [1] is employed (at next-to-leading order (NLO)). However, 
when the observable O is small and L = ln(l/(9) is large the distribution suffers 
large logarithmic contributions of the form 0{a n s L 2n ~ x ) that render the fixed order 
perturbation theory insufficient. In many cases the dominant infrared logarithms 
can be resummed using the coherent branching formalism [2, 3, 4, 5]. Furthermore, 
significant power suppressed effects are present and have been phenomenologically 
studied [6]. To accurately describe the data, all of these contributions - fixed order, 
infrared resummation and power correction - are found to be necessary. Perhaps 
the most important measurement extracted from three-jet event shape data is the 
determination of the strong coupling constant. 

Four-jet event shape observables also contain useful information about QCD. 
They are more sensitive to the triple gluon vertex and therefore the true gauge 
structure of QCD [7] and also to the presence of other light coloured particles such 
as the gluino [8] that can be pair produced by gluon splitting. However, these 
variables have received much less attention partly because they are suppressed by 
an additional power of a s requiring a second gluon to be radiated but also because 
the theoretical description is much less developed. Recently however, four separate 
general purpose Monte Carlo programs have been developed to estimate the next-to- 
leading order, 0(a 3 s ), corrections to four-jet event shapes: MENLO PARC [9], DEBRECEN 
[10] and MERCUTIO [11] employing on the one-loop helicity amplitudes for e + e~ — > 
4 partons [12] and EERAD2 [13] based on the interference of the one-loop matrix 
element with tree level [14]. For most four-jet event shape variables, the situation 
is the same as for the three-jet event shape variables; the NLO corrections are very 
large and, for renormalisation scales of the order of the centre-of-mass energy Q, still 
undershoot the data significantly [13]. This points at the presence of large infrared 
effects as well as significant power corrections. With the exception of the four-jet 
rate in certain schemes [15] and near-to-planar three-jet events [16], the issue of 
resumming infrared logarithms for specific four-jet event shape variables has not yet 
been addressed. Similarly, power suppressed effects have only been studied for the 
out-of-event plane momentum distribution [16]. It is the goal of this Letter to resum 
the leading 0[a n s L 2n ~ x ) and next-to-leading 0[a n s L 2n ~ 2 ) infrared logarithms to all 
orders in the coupling constant for the light hemisphere mass pi and narrow jet 
broadening Bn distributions. 
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To be more precise let us first define these variables. We first separate the 
event at centre-of-mass energy Q = into two hemispheres Hi, H 2 divided by the 
plane normal to the thrust axis n T . Particles that satisfy Pi-n T > are assigned 
to hemisphere Hi, while all other particles are in H 2 . Jet broadening measures the 
summed scalar momentum transverse to the thrust axis in one of the hemispheres 
while the hemisphere mass is the invariant mass of the hemisphere, 

B N = mm Pfc ' — : 1.1 

<=1 ' 2 2 EJPfc| 

1 ' E p*) ■ ( L2 ) 

PfcS-ffi / 



Pl = - ■ mm 

S i=l,2 



Note that this definition of the light hemisphere mass is the common modification of 
the original variables suggested by Clavelli [18]. These four-jet event shape variables 
are intimately connected to their three-jet event shape counterparts, the wide jet 
broadening and heavy hemisphere mass, 
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where 



that have the property of exponentiation. That is to say that the fraction of events 
where the observable O = B w or O = pn has a value less than O obeys the expo- 
nentiated form, [4, 5] 

R(0,a s (Q 2 ))= [°-^dO (1.5) 

= C(a s (Q 2 ))£(0, a s (Q 2 )) + D(0, a s (Q 2 )) (1.6) 

oo 

C(a s (g 2 )) = l + ^C n < (1.7) 

71=1 

(\ oo n+1 

S(0,a s (Q 2 ))j =^^G nm <L m 
' n=l m=l 

= Lgi(a s L) + g 2 (a s L) + a s g 3 (a s L) + ... (1.8) 

oo 

D(0,a s (g 2 )) = ^D„<. (1.9) 



?i=0 



Here C n and G„ m are constants and the perturbatively calculable coefficients D n — > 
as O — > 0. Knowledge of G nn+ i (or equivalently gi) together with G n allows 
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resummation of terms in E down to 0{a n s L 2n ~ x \ Partial inclusion of other subleading 
logarithms may be accomplished by retaining all knowledge of g2- The accuracy of 
R also depends on the number of terms known in C(a s (Q 2 ), with each known term 
giving information about two towers of logarithms. For example, C = 1 resums the 
0(a n s L 2n ) and 0(a n s L 2n ~ l ) terms in R, d resums the 0(a n s L 2n ~ 2 ) and £>«L 2n - 3 ) 
terms and so on. 

2. Coherent branching 

The coherent branching formalism allows the resummation of soft and collinear log- 
arithms due to the emission of gluons from a hard parton. As a specific example, let 
us consider the jet mass distribution J a (Q, k 2 ) as the probability of producing a final 
state jet with invariant mass k 2 from a parent parton a produced in a hard process 
at the scale Q 2 . For an initial quark this is [3] 



Y2 



1 

7& 



Ji(Q 2 ,k 2 )=5(k 2 )+ J% J dz &(z 2 (l - z) 2 q 2 - Q 2 )P™[a s (z 2 (l - z) 2 )q 2 ),z] 



J dq 2 j dk ,2 5 (k 2 - z(l - z)~q 2 - k —- J«(z 2 q 2 , k' 2 )J 9 ((l - z) 2 q 2 , q 2 ) 

(2.1) 



o 

^2 j„2 \ 



-J q (?,k 2 ) 
normalised such that 



oo 

dk 2 J a (Q 2 ,k 2 ) = 1 (2.2) 



and where the next-to-leading order q — > qg splitting kernel in the MS scheme with 
Nf flavours is 

PmK ;] = | Cf I±£( 1 + | A .) + ..., (2 .3) 

where 

K = Ca (— - — \ - 5 -N,. (2.4) 

Eq. (2.1) has a simple physical interpretation. The first term is the possibility that 
the originating parton does not emit any radiation. The transverse momentum of 
the parton is therefore unchanged. Alternatively, the quark may branch into a quark 
and a gluon subject to the phase space constraints of two body decay and which 



3 



(a) 




(b) 



Figure 1: Event pictures of (a) two-jet configurations from quark-antiquark final states 
and (b) three-jet configurations originating from qqg events where the gluon is the hardest 
parton. The cones represent coherent soft and collinear gluon emission. The thrust axis is 
denoted by a dashed line. 

subsequently undergo further emissions. This is described by the term proportional 
to J q J 9 . The last term is due to virtual corrections and ensures that soft and 
collinear singularities are regularised. A similar equation holds for J 9 and involves 
the g — > gg and g — > qq splitting kernels. Solving the integral equation for J a is 
equivalent to resumming the infrared logarithms. In particular, the probability that 
an isolated parton a forms a jet with mass less than pQ 2 is given by 



In Ref. [3], ln(E^) is solved to next-to-leading logarithmic accuracy. 

Similarly we can define the function T a (Q,\t t ;pt) [5] which describes the dis- 
tribution of the summed scalar transverse momentum p t in a jet of parton type a 
produced with vector transverse momentum k t , at scale Q. The structure is identical 
to Eq. (2.1). 

3. The probabilistic interpretation 

We can apply the coherent branching formalism to event shapes in the following 
schematic way illustrated in Fig. 1. The underlying configuration in e + e~ annihilation 
is the production of a quark-antiquark pair aligned with the thrust axis (Fig. 1(a)). 
Each parton then undergoes soft and collinear gluon emission (denoted by the open 
cone at the head of the parton). This contribution describes small angle and soft 
emission accurately in two-jet-like events when Bw and pn are small and gives rise to 
the exponentiated first term in Eq. (1.6). However, it does not describe the possibility 
of wide angle gluon emission shown in Fig. 1(b) where the gluon is the hardest 
parton. This three-jet-like matrix element correction is not logarithmically enhanced 
at 0(a s ) and first contributes at next-to-next-to-leading logarithmic accuracy and 
can therefore be neglected. 




P Q 2 



dk 2 J a (Q 2 ,k 2 ). 



(2.5) 
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Let us first focus on the hemisphere masses. To next-to-leading logarithmic 
accuracy, the fraction of events with heavy hemisphere mass less than pu i.e. k\ < 
PhQ 2 and k\ < puQ 2 is given by the two jet contribution 

R H (pH,a s (Q 2 )) = 

oo oo 

J dk\ j dk 2 2 J\Q 2 ,k\) J q (Q 2 ,k 2 2 ) Q( PH Q 2 - k\) Q(kl - k\) 



oo oo 

+ JdklJ dk 2 2 J q (Q 2 , k\) J q (Q 2 , k 2 2 ) Q(phQ 2 - kl) Q{k 2 2 - kl), (3.1) 



o o 



and where the constraint that pu 1 has suppressed the three and more jet contri- 
butions. Following the steps of Ref. [5], we can rewrite this formula using the phase 
space restrictions as 

RH( P H,a s (Q 2 )) = [j: q H (p H ,a s (Q 2 ))] 2 . (3.2) 

where 

PhQ 2 

Z« H (p H ,a s (Q 2 ))= J dk 2 J q {Q 2 , k 2 ). (3.3) 
o 

The fraction of events with light hemisphere mass less than p^ also receives 
contributions from the two-jet configuration when k\ < PlQ 2 and k\ > kl and vice 
versa. Altogether we have, 

R L (PL,a s (Q 2 )) = 



oo oo 



J dk\ j dk 2 2 J\Q 2 ,kl) J q (Q 2 ,k 2 2 ) Q(p L Q 2 - k\) Q(k 2 2 - k\) 



oo oo 

+ JdklJ dk 2 2 J%Q 2 , kl) Ji(Q 2 , k 2 2 ) Q(plQ 2 - kl) Q(kl - k\). (3.4) 
o o 

Simplifying the phase space constraints, we find 

R L (PL,a s (Q 2 )) = 2^ H (p L ,a s (Q 2 ))- [£^(p L , a s (Q 2 ))] 2 . (3.5) 

The functions that resum the logarithms for the light hemisphere mass are the same 
as those that resum the logarithms for the heavy hemisphere mass. Now however, 
exponentiation in its purest form is spoiled because the final result is a sum of terms. 
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The analysis for the narrow jet broadening proceeds in the same way. We find 
that Rn, the probability of finding an event with a light hemisphere mass less than 
-Bat, is given by 



R N (B N ,a s (Q 2 )) 



-Biv<l 



2 Zl(B N ,a s (Q 2 )) - [Zl(B N ,a s (Q 2 ))} 2 (3.6) 



where the probability of obtaining a jet with summed scalar transverse momentum 
p t with respect to the jet axis less than 2BQ starting from a parton of type a, S^, 
is given by 



4. All-orders resummation of large logarithms 

In this section we discuss the all-orders resummation of leading logarithms 0{a n s L 2n ) 
and next-to-leading 0{a n s L 2n ~ l ) logarithms to all orders in the coupling constant. 
From eqs. (3.5) and (3.6), we see that to determine Rl and R N requires knowledge 
of Y> a H and respectively. Both of these functions have the exponentiated form of 
Eq. (1.6) and the corresponding functions gi(a s L) and g2(a s L) determined. 

Explicit expressions for valid to this order are given in [3] and, introducing the 
renormalisation scale dependence in the standard manner and dropping the parton 
index, we reproduce them here for illustrative purposes, 



(-2BQ 



w 



(B,a s (Q 2 )= / T a (Q,0,p t )dp, 



't- 



(3.7) 



Jo 



£tf(p,a s (^ 2 ),^) 



exp[^(q s (/i 2 ),L)] 
T[l-S(a s (^),L)] 

expfL/^x) + h{x) + x 2 f[{x) Hfi 2 /Q 2 )] 



+ OWL 



)(4-l) 




where 



L = ln(l/p), 



x = (3 a s (fi 2 )L. 



(4.2) 
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The functions f\, f 2 and f[ are 



f2(x) 



A W 



2irf3ox 
AW 



(1 - 2x) ln(l - 2x) - 2(1 - x) ln(l - x) 
2m(l-x) -m(l-2x 



(4.3) 



7T/? 



2nf3$ 

AW 
2irf3 x 2 



ln(l - x) - ln(l - 2x) 
1 



ln(l - 2x) - 2 ln(l - x) + - ln 2 (l - 2x) - ln 2 (l - x) 

2 



ln(l - 2x) - 21n(l - x) 



(4.4) 
(4.5) 



with 



For quarks, 



A) 



11C A -2A^ 
12tt ! 



17CU 2 - 5C A iV/ - 3C F jV/ 
24^2 • 



= C f 



aw = IcfK, BW = - 3 - Cf , 



(4.6) 



(4.7) 



with K given by Eq. (2.4). 

Altogether Eqs. (4.1) to (4.7) are sufficient to determine Rn of Eq. (3.6) to next- 
to- leading logarithmic accuracy. The expression for therefore correctly sums to 
all orders in the strong coupling, a s , only the leading two towers of large logarithms 
from 0(a n s L 2n ) down to 0(a n s L 2n - 1 ). 

The resummed formula (3.5) does not include the 0(a 2 L 2 ) or 0(a 2 L) terms 
(just as the analogous formulae for resumming three jet variables do not include 
the 0(a 2 s L 2 ) or 0(a 2 L) terms) present in the lowest order perturbative coefficient. 
Similarly, it does not produce the 0(alL 4 ) to 0(a^L) terms that occur in the next- 
to- leading order perturbative coefficient. The perturbative calculation provides the 
a 2 and contributions exactly and therefore the most significant omitted term is 
0(a 4 s L 6 ), see Fig. 2. 

Precisely the same discussion applies to the narrow jet broadening. Using the 
coherent branching formalism the functions gi(a s L) and g 2 (a s L) have been deter- 
mined for Ti W and Catani et al (CTW) have provided analogous expressions for 
T, w that are given in [5]. However, in doing so certain simplifying approximations 
concerning the recoil transverse momentum have been made. Dokshitzer and collab- 
orators [19] have found that treating the quark recoil more carefully causes the CTW 
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Figure 2: The towers of infrared logarithms appearing in the four-jet event shape rates Rl 
and Rn- The resummation includes the leading and next-to-leading logarithms denoted 
by filled squares and the complete 0{a 2 s ,o? s ) contributions from fixed order denoted by the 
empty squares. All other terms are incomplete and denoted by empty circles. The black 
filled squares denote terms generated purely in the two-jet limit. 



result for T. w to be adjusted by a multiplicative factor. This modifies g2(a s L) but 
leaves Gu unchanged. The same correction is in principle needed for the narrow jet 
broadening but is relevant beyond the next-to-leading logarithmic accuracy. Never- 
theless, we chose to apply the correction to the CTW result. We do not give the final 
form for here, but instead refer the interested reader to Ref. [19]. Inserting the 
recoil-corrected form for in the resummed expression for (3.6) again allows 
resummation of the leading and next-to- leading logarithms as shown in Fig. 2. 

5. Numerical results 

As usual, the resummed result contains part of the fixed order perturbative contribu- 
tion and the overlap must be removed by matching. This is done by expanding the 
resummed result as a series in the strong coupling constant and explicitly removing 
the terms corresponding to the fixed order calculations. This can be achieved in 
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several ways of which R matching and ln(i?) matching are the most common. In 
the R matching scheme, the coefficients of each of the unsummed logarithms present 
in the fixed order perturbative coefficients must be numerically extracted. For the 
four-jet event shape observables discussed here, this corresponds to determining the 
coefficient of a 2 L from the lowest order perturbative contribution and the coefficients 
of a^L 2 and a z s L from the next-to-leading order contribution. This is impractical. 
However, in the more commonly used ln(R) matching scheme, it is assumed that 
the fixed order result exponentiates and therefore it is not necessary to make this 
extraction because any logarithmic terms remaining after subtracting the overlap 
from the fixed order contribution are exponentially suppressed. Of course, for the 
four-jet variables considered here, R does not exponentiate. Nevertheless, in the 
ln(i?) matching procedure, any remaining logarithmic terms are still exponentially 
suppressed. We therefore employ the ln(R) matching procedure. 
To be more precise, the resummed result has the form 

Rresummed = Ci2S — C 2 £ 2 + 0(a s ), (5.1) 

where Cj are power series in a s and 0(a s ) represents the contribution from three- 
jet configurations. The 0(a s ) term of C\ and C 2 are related by the constraint that 
there can be no logarithmic contributions at 0(a s ) - the expansion of R containing 
logarithms starts at 0(a 2 ). However, to the next-to-leading logarithmic accuracy 
relevant for this paper, we can set C± — C 2 — 1. The InR matching formula then 
dictates that 

ln(R) = hl{Rresummed) + ^{R fixed) ~ ^(Rexpanded) , (5-2) 

where In Rf ixec i and \nR expan d e d are obtained by expanding the fixed order result and 
the In Rresummed through to 0{ofy respectively. Using the expansion of ln(E) as a 
series, 

oo n+1 

= E E Gnm^L™, 
n=l m=l 

and with C\ — C 2 — 1, we see that, 

^ n (Rexpanded) = —{GnL + Gi2L 2 ) 2 a 2 

-2(G U L + G 12 L 2 )(G 2 iL + G 22 L 2 + G 23 L 3 + i(G n L + G 12 L 2 ) 2 )a 3 s 
+0(at). (5.3) 

At next-to-leading logarithmic accuracy, logarithmic terms from the fixed order con- 
tribution remain in ln(i2) given by Eq. (5.2), however, they are exponentially sup- 
pressed when translated back to R. 

We expect that at large values of the observable O the resummed result is dom- 
inated by the fixed order calculation. However, the resummed expressions (3.5) and 
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(3.6) valid in the small pl and limits do not contain information about the kine- 
matic endpoints of the distributions. To ensure that the resummed result vanishes 
at the endpoint we make the substitution 

h^h~o^ +1 - (5 - 4) 

where O max corresponds to the endpoint of the distribution at the accuracy of the 
fixed order calculation. We use p™ ax = 0.167 and 5™ ax = 0.204. 

Numerical results for the light hemisphere mass and for the narrow jet broadening 
are shown in figures 3 and 4 respectively. Throughout we set p = Q = M z and use 
a s (M z ) = 0.118 corresponding to the current world average. The next-to-leading 
order result which diverges at small values of the event shapes is taken from [13] and 
is evaluated at p = Q. 

Figures 3 and 4 show that the resummations are extremely important for pi < 
0.01 and B^ < 0.02. Rather than the divergent fixed order prediction, we have 
the more physical resummed result that the probability of finding events with no 
radiation (very small values of pl and B N ) are vanishingly small. For the reference 
value of p — Q, the peak position of the pi distribution occurs at pi = 0.01 with 
a height of 0.33, while for B N the peak occurs at B N = 0.02 with a value of 0.48. 
At larger values, the resummation changes the NLO prediction by a more moderate 
amount indicating that uncalculated higher order corrections are under control. At 
very large values of O, the resummed and NLO predictions coincide because of the 
matching procedure of Eq. (5.4). 

To illustrate the residual renormalisation scale dependence, we also show the 
effect of varying p by a factor of 2 either side of the reference value p = Q. We see 
that around the peak region, different scale choices alter the prediction by ±10%. 

The infrared resummation significantly improves the perturbative prediction for 
the event shape observable. However, when comparing with experimental data 
we should be aware that important non-perturbative hadronisation corrections are 
present. The effect of hadronisation on the distribution is to shift the value of the 
observable away from the two-jet region, 

0^0 + NP , (5.5) 

where the non-perturbative correction depends on the typical hadron scale 0(1 GeV) 
and is suppressed by a power of Q. In principle these power corrections can be esti- 
mated using the dispersive approach of Ref. [20] where a non-perturbative parameter 
pi is introduced to describe the running of a s in the infrared region. For the asso- 
ciated three jet variables the non-perturbative corrections are typically estimated 
to be 0(1 GeV/Q) for p H [20] and O(0.3 GeVln(l /B w )/Q) for B w [21] and arise 
through the hadronisation of one of the two jets in the event. Because the four-jet 
event shapes are largely related to what happens in the second jet, we might expect 
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Figure 3: The resummed (solid) and fixed order NLO (dashed) predictions for the light 
hemisphere mass distribution at \i = Q = Mz- The resummed prediction at // = Q/2 

(fi = 2Q) is shown as a dotted (dot-dashed) curve. 



n 1 i 1 i 1 i 1 i 1 i 1 i 1 i 1 i r 

.*'"■•. Resummed inc. recoil at u.=Q/2 




Figure 4: The resummed (solid) and fixed order NLO (dashed) predictions for the narrow 
jet broadening distribution at fx = Q = Mz- The resummed prediction at fi = Q/2 

(fi = 2Q) is shown as a dotted (dot-dashed) curve. 



that the hadronisation corrections are similar. To illustrate the potential effects of 
hadronisation in the four-jet event shapes, we just transfer these corrections directly 
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Figure 5: The resummed prediction for the light hemisphere mass distribution a t 
fi = Q = Mz modified by a non-perturbative power correction pi — ► pi + 1 GeV/Q. The 
resummed prediction without power correction is shown as a dashed line. For compari- 
son, we also show the charged hadron data collected at the Z resonance by the DELPHI 
collaboration [22]. 



so that, 



Pl -> Pl + 
B N — »■ B N + 



1 GeV 



Q ' 

0.3 GeVln(l/ J B jV ) 

Q 



(5.6) 
(5.7) 



The simplified hadronisation correction applied to the resummed distributions for the 
light hemisphere mass and narrow jet broadening with /i — Q is shown in figures 5 
and 6 respectively. To emphasize the dramatic effect the power correction has, we 
also show the uncorrected predictions. There are two effects. First the distribution 
is shifted to the right by an amount Onp and second, the distribution is rescaled by 
a factor (O + Onp) / O. In the region of the turnover where O is of the same order as 
Onp there is an enhancement of (9(100%). The hadronisation correction is smaller 
at larger values of O. 

For comparison, we also show the charged hadron data collected by the DELPHI 
Collaboration [22] at the Z resonance. We see remarkable agreement (strikingly so in 
view of the simplified hadronisation correction applied here). The only discrepancy 
is at very small values of O < Onp where individual hadrons in the light/narrow 
hemisphere will significantly affect the value of O. We do not expect to successfully 
describe such events. 
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Figure 6: The resummed prediction for the narrow jet broadening distribution —jf--^§^ 
at fi = Q = Mz modified by a non-perturbative power correction Bn — > + 
0.3 GeVm(l/_Bjv)/Q- The resummed prediction without power correction is shown as 
a dashed line. For comparison, we also show the charged hadron data collected at the Z 
resonance by the DELPHI collaboration [22]. 



6. Conclusions 



In this paper, we presented predictions for the light hemisphere mass and the nar- 
row jet broadening distributions where the infrared logarithms have been resummed 
to next-to-leading logarithmic order using the coherent branching formalism. The 
resummed expressions do not exponentiate, but involve the difference of exponential 
factors so that the perturbative series starts at 0(a 2 ). The expressions presented 
in sections 3 and 4 resum the leading 0{a n s L 2n ' v ) and next-to-leading 0{a n s L 2n ~ 2 ) 
infrared logarithms to all orders in the coupling constant. Taken in conjunction with 
the fixed order 0(a 3 s ) non- logarithmic terms, the first neglected term is of 0(a i s L % ). 
The numerical results shown in Figs. 3 and 4 indicate that the resummation effects 
are sizeable at small values of the event shape parameter and that the resummation 
procedure significantly improves the perturbative prediction. 

However, because important non-perturbative hadronisation corrections are 
present, resummation alone is insufficient to describe the experimental data. Ap- 
plying simple power corrections similar to those obtained for the heavy hemisphere 
mass and the wide jet broadening gives good qualitative agreement with the avail- 
able data from LEP. We anticipate that the improved theoretical description of the 
four-jet event shape distributions presented here can be combined with a more so- 
phisticated hadronisation correction based on the dispersive approach of Ref. [20]. 
The data from LEP can then be used to further test the structure of QCD in four-jet 
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like events and extract values of the strong coupling constant. 
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